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' Abstract 

' In this work we study a generalized nonlocal thermistor problem with fractional-order 

Riemann-Liouville derivative. Making use of fixed-point theory, we obtain existence and 
uniqueness of a positive solution. 
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1 Introduction 

Joule heating is generated by the resistance of materials to electrical current and is present in 
any electrical conductor operating at normal temperatures. The heating of such conductors has 
undesirable side effects. Problems dealing with the combined heat and current flows were consid- 
ered in [511^ 1 lllfTB] , where various aspects of the so-called thermistor problem were analyzed. The 
mathematical model of the nonlocal steady thermistor problem has the form 

Au=^i^, (1) 

where A is the Laplacian with respect to the spacial variables. Such problems arise in many 
applications, for instance, in studying the heat transfer in a resistor device whose electrical con- 
ductivity / is strongly dependent on the temperature u. The equation ([1]) describes the diffusion 
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of the temperature with the presence of a nonlocal term as a result of Joule effect. Constant A 
is a dimensionless parameter, which can be identified with the square of the applied potential 
difference at the ends of the conductor. Function u represents the temperature generated by the 
electric current fiowing through a conductor. For more description, we refer to |12[|21j . A deep 
discussion about the history of thermistors, and more detailed accounts of their advantages and 
applications to industry, can be found in jTll[13]. In [4] Antontsev and Chipot studied existence 
and regularity of weak solutions to the thermistor problem under the condition that the electrical 
conductivity f{u) is bounded. 

Fractional differential equations are a generalization of ordinary differential equations and 
integration to arbitrary noninteger orders. The origin of fractional calculus goes back to Newton 
and Leibniz in the seventieth century. In recent years there has been a great deal of interest in 
fractional differential equations. They provide a powerful tool for modeling and solving various 
problems in various fields: physics, mechanics, engineering, electrochemistry, economics, visco- 
elasticity, feedback amplifiers, electrical circuits and fractional multipoles — see, for example, [HI 
|8|10, 18-20 and references therein. Using fixed point theorems, like Shauder's fixed point theorem, 
and Banach's contraction mapping principle, many results of existence have been obtained to linear 
and nonlinear equations and, more recently, to fractional derivative equations. The interested 
reader can see [2l[3]. For a physical meaning to the initial conditions of fractional differential 
equations with Riemann-Liouville derivatives we refer to [Tl lTClfTB] . 

Our main concern in this paper is to prove existence and uniqueness of solution to a general 
fractional order nonlocal thermistor problem of the form 

D^'-u^^f^^ + h{t), te{0,T), 

U!f{u)dxf (2) 

/^u(i)|t=o = 0, V/?e(0,l], 

under suitable conditions on / and h (see Theorem 13. 2p . We assume that T is a fixed positive 
real and a > a parameter describing the order of the fractional derivative. In the literature 
we may find a great number of definitions of fractional derivatives. In this paper, the fractional 
derivative is considered in the Riemann-Liouville sense. In the case a = 1 and h — 0, the fractional 
equation ([2|) becomes the one-dimensional nonlocal steady state thermistor problem. The values 
of < a < i correspond to intermediate processes. We further prove the boundedness of u (see 
Theorem 13. 3[) . which is of considerable importance from a practical and physical point of view: it 
is interesting to keep the temperature from exceeding some extremal values that may damage the 
conductor. 



2 Preliminaries 

In this section, we give some basic definitions and preliminary facts that are used further in the 
paper. Let < a < ^ and X = (C([0, T]), || • ||), where C([0,r]) is the space of all continuous 
functions on [0,T]. For x G C([0,T]), define the norm 

\\x\\ = sup {e'^'^xm, 
te[o,T] 

which is equivalent to the standard supremum norm for / G C{[0,T]). It is used in literature in 
many papers, see for example [1 . The use of this norm is technical and allow us to simplify the 
integral calculus. By L^{[0,T],M.) we denote the set of Lebesgue integrable functions on [0,T]. 
Throughout the text c denote constants which may change at each occurrence. As in ,4,, we 
consider that the electrical conductivity is bounded. We now assume the following assumptions: 

(HI) / : M+ IR+ is a Lipschitz continuous function with Lipschitz constant Lf such that 
ci < /(li) ^ C2, with ci,C2 two positive constants; 

(H2) h is continuous on (0,T) with h e L°^{Q,T). 
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Definition 2.1. The fractional (arbitrary) integral of order a G IR+ of a function / G Li[a,b] is 
defined by 



where F is the gamma function (see, e.g., [14l[T6l[T7lfT9] ) . For a = we put /" :— Ig. 

Remark 2.1. For /, g G Li[a,b] one has 

Note also that I°'f{t) G C(R+) for / G C(M+) and, moreover, /"/(O) = 0. 

Definition 2.2 (see, e.g., P^[TB1[T71[T5] ) . The Riemann-Liouville fractional (arbitrary) derivative 
of order a G (n — 1, n), n G N, of function / is defined by 

3 Main Results 

Our main result asserts existence of a unique solution to on C(M+) of the form 

^ ' ' ^ ^ + h{s) } ds . 



(3) 



F(2a) [{J^f{u)dxy 



3.1 Existence and Uniqueness 

We begin by proving the equivalence between ([2]) and ([3]) on the space C(R+). This restriction 
of the space of functions allows to exclude from the proof a stationary function with Riemann- 
Liouville derivative of order 2a equal to d ■ t^"'^, d G M, which belongs to the space Ci_2Q[0,r] 
of continuous weighted functions. 

Lemma 3.1. Suppose that a G (0, ^). Then the nonlocal problem ^ is equivalent to the integral 
equation ^ on the space C(W^). 

Proof. First we prove that ^ implies For t > equation ([5]) can be written as 

^^^"^ , +Ht). 

'^^ {lof{u)dx] 
Integrating both sides of the above equation, we obtain 



l'-"-u{t) - l'-"^u{t)\t^o = / < ' ,2 + Hs) } ds. 

^0 [(j^f{u)dx) 



Since < 1 - 2q; < 1, 



I^-^°^u{t)^ { ^h{s))ds. 



' i(/o'/(") 



dx 
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Applying the operator / to both sides, we get 



Iu{t) = / 



Differentiating both sides, 



Jo f{u)dx 
A/(«) 



+ h{t) 



hit) 



Let us now prove that ([3|) iniphes ([2]). Since u e C{I) and G C(/), applying the operator 

to both sides of ((21) one obtains 



Jq fiu)dx 



hit) 



= I 



A/(w) 



Jo f{u) dx 



hit) 



Differentiating both sides of the above equality. 



Dl'^-^'^uit) = DI 



A/(u) 



Jo /(") 



hit) 



Then, 



A/(ii) 



(/o^ /(") 



hit), t>0. 



□ 



Theorem 3.2. Let f and h satisfy hypotheses (-ffl) and iH2). Then there exists a unique solution 
u e X of ^ for alio < X< — . 



Proof Let F : X ^ X he defined by 



Fu = r 



A/H 



Then, 

\Fu - Fv\ = 



j2a I 



A/(zi) 



Jo f{u)dx 



hit) 



[j^fiu)dx) {jofiv)dx) 



< 



T2a 



r2a 



[jo /(") '^^ 
A 

{jo /(") '^^ 



Afiu)-fiv)) + Xfiv) 



[[j^fiu)dx) {jofiv)dxy 



-ifiu)-fiv)) 



{jofiu)dx) {jofiv)dx) 

(4) 



4 



We estimate each term on the right hand side of (|1]) separately. Using then the fact that / is 
Lipshitzian, we have 



r2a 



[lo /(") dx 



-if{u)-f{v)) 



< 



1 



< 



(ciT)2 

1 
1 



Xl'"{\f{u)-fiv)\} 
XLfl''^{\u-v\} 



(5) 



(ciT) 



{t~s) 



2q-1 



r(2a) 



\u{s) — v{s)\ds. 



Since 



/o r(2a) 
it follows from (O that 



Nt 2a-l 1 

e-''dr < 



-Nt 



I' 



r(2a) Jo 
-{f{u)-f{v)) 



r(2a) 



1 



2a-l 



[ciTy 7o r(2a) 



< 



< 



Jo r(2Q!) 

AL/ sup{e-^*|u(i) - v{t)\} r ^^^^^e-^(*-^)ds 



(cir)2 
1 



< 



< 



1 



1 



r(2a) 



r(2«) 



2^ |-Aft ^2q-1 



(cir)2 "iV2"yo r(2a) 

On the other hand, similar arguments as above yield to 



1 



1 



r2a 



A/(«) 



f{u) dx 



f{v) dx 



C2 , 



< 



(cir)4 

2c?T 



r2a 



f{u)dx\ - f{v)dx 



{f{u)-I{v))dx\ / + 



- (cir)4 ^ 



— v{x)\ dx 



- Mr^'^^X r(2a) 



— ^(a;)! do: ds. 



(6) 
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Then, 




^(/o fHdxy (^J^ f{v)dx 




{t-sf 



r(2a) 



e^^(*~^)e^^^|M(x) - v{x)\dx ds 



< 



< 



< 



< 



< 



< 



< 



2clT 

Wry 

2c|T 

2cIT 
(ciT) 



XLf sup{e~'^>(i) -^(01} 



XLf sup{e-^'\u{t)-v{t)\} 



r(2a) ' 



-Aft 



e-^(*-") dx I ds 
1 



XLf\\u~v\\ 



2q-1 



-Nt 



N 

„NT\t 



N 



N 



7V2Q+1 



U — V 



U — V 



r(2a) 



r(2a) 



1) ds 



(7) 



\2a-l 



r(2a) 
r(2a) 



r(2a) 



Gathering ©-Q, we get 

e-^^\Fu-Fv\ < 

Then, we have 



1 



(cry 



\\Fu - Fv\\ < 



( + (ciT)4 



^2q 



Choosing A > such that (yjj^ 



TP 



(ciT)4 



< 1, the map F : X X is Si contraction 



and it has a fixed point u = Fu. Hence, there exists a unique u E X that is the solution to the 
integral equation ([3]). The result follows from Lemma [3TT] □ 

3.2 Boundedness 

We now show that the condition that the electrical conductivity f{u) is bounded (hypothesis (HI)) 
allows to assert boundedness of u. 

Theorem 3.3. Under hypotheses {HI) and {H2) and A > 0, if u is the solution of then 

.. .. ((^/(0)+/lo 



u < 
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Proof. One has 



{cT^Jo T{2a) ' V(ciT)2-'^ ' ' "^JJo r(2a) 

Then, 

rt u „^2a-l 



-Art 



-{c^ryk r(2a) msjids + ^^^^^^^/W + ftoojy^ j,^2a) ^ ^ ''^ 

(ciT)^ Jo r(2a) [ciTy Jo r(2a) 

-(cir)2X r(2a) 1"^'^!^'+ 7V2„ r(2a)' 

Using Gronwall's lemma, we have 

< '■^ p Y^ia Jo r(2c«) 

^ ^P2^ 

Then, 

„ „ (7sitf/(o) + /loo) (tstW/W + /ioo) ^ 

" " — ]\l2a ]\f2a 

and we conclude that u is bounded. □ 



Acknowledgements 

Work supported by FEDER funds through COMPETE — Operational Programme Factors of 

Competitiveness ( "Programa Opcracional Factorcs de Compctitividadc" ) and by Portuguese funds 
through the Center for Research and Development in Mathematics and Applications (University 
of Aveiro) and the Portuguese Foundation for Science and Technology ("FCT — Fundagao para 
a Ciencia e a Tecnologia" ) , within project PEst-C/MAT/UI4106/2011 with COMPETE number 
FCOMP-01-0124-FEDER-022690. The authors were also supported by the project New Explo- 
rations in Control Theory Through Advanced Research (NECTAR) cofinanced by FCT, Portugal, 
and the Centre National de la Recherche Scientifique et Technique (CNRST), Morocco. 

The authors are grateful to an anonymous referee for several corrections and useful comments. 



7 



References 

[1] S. Abbas, Existence of solutions to fractional order ordinary and delay differential equations and 
applications, Electron. J. Differential Equations 2011, No. 9, If pp. 

[2] R. P. Agarwal, M. Benchohra and S. Hamani, Boundary value problems for fractional differential 
equations, Georgian Math. J. 16 (2009), no. 3, 401-4ff. 

[3] R. P. Agarwal, Y. Zhou and Y. He, Existence of fractional neutral functional differential equations, 
Comput. Math. Appl. 59 (20f0), no. 3, f095-ff00. 

[4] S. N. Antontsev and M. Chipot, The thermistor problem: existence, smoothness uniqueness, blowup, 
SIAM J. Math. Anal. 25 (f994), no. 4, f f28-ff56. 

[5] G. Cimatti, Existence of weak solutions for the nonstationary problem of the joule heating of a 
conductor, Ann. Mat. Pura Appl. (4) 162 (1992), 33-42. 

[6] L. Gaul, P. Klein and S. Kempfle, Damping description involving fractional operators, Mech. Syst. 
Signal Process. 5 (1991), 8f-88. 

[7] C. Giannantoni, The problem of the initial conditions and their physical meaning in linear differential 
equations of fractional order, Appl. Math. Comput. 141 (2003), no. 1, 87-f02. 

[8] R. Hilfer, Applications of fractional calculus in physics, World Sci. Publishing, River Edge, NJ, 2000. 

[9] S. D. Howison, J. F. Rodrigues and M. Shillor, Stationary solutions to the thermistor problem, J. 
Math. Anal. Appl. 174 (1993), no. 2, 573-588. 

[f 0] A. A. Kilbas, H. M. Srivastava and J. J. Trujillo, Theory and applications of fractional differential 
equations, North-ffoUand Mathematics Studies, 204, Elsevier, Amsterdam, 2006. 

[f f] K. Kwok, Complete guide to semiconductor devices, McGraw-Hill, New york, 1995. 

[12] A. A. Lacey, Thermal runaway in a non-local problem modelling Ohmic heating. H. General proof of 
blow-up and asymptotics of runaway, European J. Appl. Math. 6 (1995), no. 3, 201-224. 

[13] E. D. Maclen, Thermistors, Electrochemical publication, Glasgow, 1979. 

[14] K. S. Miller and B. Ross, An introduction to the fractional calculus and fractional differential equa- 
tions, A Wiley-Interscience Publication, Wiley, New York, 1993. 

[15] D. Mozyrska and D. F. M. Torres, Modified optimal energy and initial memory of fractional 
continuous-time linear systems. Signal Process. 91 (2011), no. 3, 379-385. lairXiv : 1007 ■ 3946] 

[16] 1. Podlubny, Fractional differential equations. Mathematics in Science and Engineering, 198, Academic 
Press, San Diego, CA, 1999. 

[17] 1. Podlubny and A. M. A. El-Sayed, On two definitions of fractional calculus, preprint UEF (ISBN 
80-7099-252-2), Solvak Academy of Science-Institute of Experimental Phys (1996), 03-06. 

[18] J. Sabatier, O. P. Agrawal, J. A. Tenreiro Machado, Advances in fractional calculus. Springer, Dor- 
drecht, 2007. 

[19] S. G. Samko, A. A. Kilbas and O. I. Marichev, Fractional integrals and derivatives, translated from 
the 1987 Russian original, Gordon and Breach, Yverdon, 1993. 

[20] H. M. Srivastava and R. K. Saxena, Operators of fractional integration and their applications, Appl. 
Math. Comput. 118 (2001), no. 1, 1-52. 

[21] D. E. Tzanetis, Blow-up of radially symmetric solutions of a non-local problem modelling Ohmic 
heating. Electron. J. Differential Equations 2002 (2002), no. 11, 26 pp. 



8 



